How low is the thermodynamical limit? 
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We analyze how can some dynamical process lead to exponential (Boltzmannian) distribution 
without instantaneous thermal equilibrium with a heat bath. We present a model for parton dressing 
which re-combines the exponential from cut power law minijet distributions fitted to pQCD as a 
-j: • limiting distribution. Thermal models and power law tails are interpreted within the framework of 

\ the more general Tsallis distribution. 
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O ■ I. INTRODUCTION 

:z;: 

' Not only in heavy ion collisions, but also in pp and even in e^e collisions exponential transverse mass spectra has 
^s^j ' been observed besides a power law tail at relatively high transverse momenta. The exponential shape at medium pT 
was often interpreted as an indication that the source of the observed hadrons can be described by a temperature 
' related to the slope of these spectra. This view is most prominent in the thermal modelsQ which assume a thermally 
K*" , equilibrated state, well described by the familiar equilibrium thermodynamics. 

' The experimental data from CERN SPS and RHIC show an exponential shape in the transverse mass spectrum 
, best fitted by exp((m — mT)/T), where rriT = \J vn? + pi^ is the transverse mass of a particle with mass m^]. The 
' stiffness of these spectra, T, has still to be corrected for collective flow effects: for high momenta a relativistic 
^-H ! blue-shift factor, for low momenta a mass dependent kinetic flow energy addition occurs in this parameter. The 
conjectured temperature, Tq, of a thermally equilibrated state is therefore usually higher than the directly observed 
' slope parameter, T . Characteristic values at RHIC are T « 140 MeV and Tq « 210 MeV, so the assumed temperature 
(— I is well above the color deconfinement temperature predicted by recent lattice QCD calculations. 
O I There is, however, doubt whether a thermally equilibrated state can be formed in such a short time, during 
' which the final state fireball at these high bombarding energies develops. There is no known hadronic process 
which would lead to a thermal state from a far off-equilibrium initial situation on a time scale of 1 — 2 fm/c. It is 
(— I ^ strongly tempting to consider quark level mechanisms as reason for the observed exponential spectra. Several such 
• • < mechanisms has already been suggested on a subhadronic level. Ref.|3| considered the Schwinger mechanism with a 
. !^ Gaussian fluctuating string tension, which leads to exponential rriT spectra. Having color ropes in mind, the effective 
. string tension may have a Gaussian distribution due to a random color chargation process in the limit of high color. 
' Another idea is to consider the fluctuation of temperature, distributed according to a Gamma distribution In this 
5^ , case the average of exponentials leads to a Tsallis distribution, which is still well approximated by an exponential 
at soft momenta, while ends in a power-law tail at high momenta. In this paper we shall have a closer look on a 
recently suggested hadronization mechanism, on the parton recombination , as a possible source of the experimentally 
observed exponential spectra. In this case we begin with cut power-law minijet distributions Q, havinj^ a form of the 
canonical Tsallis distribution and combine their limiting distribution for many- fold recombination 8]. 

According to this goal first we review how the Gaussian distribution emerges as a limiting distribution and recall the 
central limit theorem of statistics. Then we consider a simple example not falling under the conditions of this renown 
theorem, but still leading to a non-trivial limiting distribution. Finally we relate this result to parton recombination 
and to the canonical Tsallis distribution, which gives an improved quantitative description of experimental findings 
compared to the naive thermal model. 
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II. THE CENTRAL LIMIT THEOREM 



In order to demonstrate, how fast a Gaussian limiting distribution may develop, let us first consider the distribution 
of the scaled sum of imiform random deviates. Let Xi be a uniform random deviate in (—1, 1) and Pn{x) the distribution 

of a; = \ — yi" , Xi. The statement of the cited theorem is that in the limit of infinite n the distribution of the 
variable x exactly becomes the normal distribution. 

P^{x) = e--"l\ (1) 

Before we review how and why this happens it is educating to make some guesses how large n actually should be 
in order to mistake such a distribution for a Gaussian. The answer is surprisingly low (cf. Fig^. The n = 1 case 
is the uniform distribution, the n = 2 has a triangular shape. The finiteness of the sample (m — 200000) smoothes 
the differences out, it is therefore not easy to distinguish already the n = 3 or n = 4 case from the Gaussian in a 
numerical experiment. 

Now we turn to the proof of the limiting distribution being a Gaussian. Consider the Fourier transform of 

P„(fc) = fdx^"^^ n (1 Cdx\ E{x^S^-Y^x,). (2) 
After X- integration it turns to be a product of n equal factors: 



= ( ""^'fZ^M • (3) 
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The large n limit of this expression is a Gaussian: 

P^(fc) = hm 1-^ + ... ^ e-^'l\ (4) 

n— XX) y In J 

Fourier transforming back to x-space gives the required statement. The phenomenon is more general. Let a; be a 
weighted sum of n independent random variables: 

/n n 
(Y[w{xi)dxi) S{x- a„^Xj). (5) 
i=i j=i 

each distributed according to the same probability density w{x). Then the so called central moments, which describe 
correlations, 



(6) 



fe=0 



can be obtained easily from the central moments Cj(l) of the Pi{x) = distributions using logP„(fc) = 

nlogiZ;(fca„): Cj{n) = na^Cj(l). Whenever a„ cx Cj{n) oc n}^^/'^Cj{l) — > as n ^ oo for j > 2. This 

completes the proff of the Central Limit Theorem. 

There are exemptions to this theorem. Whenever log Poo (fc) is not differentiable arbirary times in the vicinity of 
fc = 0, or not the proper scaling factor a„ is used, or the base distribution w(x) has a long tail not approaching zero 
fast enough. In the last case there can be a limiting distribution, it is just not Gaussian. As an interesting example 
we study the convolution of normalized Lorentzians 

- 7 x^^ 

The logarithm of the Fourier transform is not differentiable at /c = 0: log ?«(/;) = — |fc|P. The n-fold convolution's 
Fourier transform also has a non-diffcrcntiable logarithm, logP„(A;) = — nP|fc|a„. It can, however, be found a finite 
limiting distribution for (n — > oo) if a„ = 1/n. In this case the Lorentzian is a fixed point of the convolution and 
scaling. The limiting Lorentzian distribution is also long-tailed, so the central moments diverge. Aiming at physical 
applications we have to regularize. We use a small parameter, m, for this purpose and at the end take the nontrivial 
double limit m ^ and n ^ oo. Starting with regularized Lorentzians, the Fourier transform has the form of the 
experimentally observed m — rriT exponential. 



w{k) = exp( ) (8) 

with T = Now logw(fc) = R{m— y/W+m^) with vanishing odd central moments, and with regular even central 
moments, C2j(l) cx ^^-Ij-i ■ The limiting distribution is obtained by keeping M = nm constant and scaling with 
a„ = 1/n again. It has the following Fourier spectrum 



P^(fc) = e^M-VW+m)/T_ 

This resembles the original regularized distribution, but with an eventually finite mass M. The low and high momen- 
tum limits are Gaussian and pure exponential, respectively. 



III. PARTON DRESSING TO EXPONENTIAL SPECTRA 



Combining hadrons from many partons also results in a limiting distribution, which is not Gaussian. We consider 
an additive quark model where the hadron mass is made of n approximatevly equal partons: M = J2^i = J^"^- The 

variable, x = ^~f'^' = -JZxi is the center of energy coordinate. Its limiting distribution, P^o (x) can be viewed as 

the hadron wave function squared. Its Fourier transform, P{k) is proportional to the hadron spectrum. 

While constructing from "equal-right" partons, we have Pn{k) = w(fc/n)". The particular cases n = 2 and 
n = 3 of this formula have been used in parton recombination models of hadronization recently. The partons before 
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hadronization show a cut power law transverse momentum distribution, as fitted to minijets obtained in the framework 
of pQCD based calculations: 

= a (l + + mV&) (10) 

Here h and c are parameters of the parton distribution, in no way connected to a thermal state, m is the transverse 
mass at zero rapidity, and finally a can be obtained from the normalization condition, wifi) = 1. Note that this 
distribution is a power law for large transverse momenta lim|;j|-^^ w{k) = ab'^ \k\~^- 

Wc construct a constituent quark and eventually hadrons from n such partons, keeping the total mass M = nm 
fixed. We arrive at 

P„(fc) = a"fl + — x/fc2 + M2j (11) 
In the n oo limit this gives the exponential distribution in the transverse mass 

P^{k) = exp ( j . (12) 

This is exponential of the reduced transverse mass Mt — M = + — M has a slope T = b/c. This is 
not a temperature, it is not related to any energy exchange with a heat reservoir. The large parameter by which 
one approaches the limiting distribution is the foldness of the recombination; in the present model the number of 
independent primordial partons dressing to a hadron constituent. 



IV. TSALLIS DISTRIBUTION AND THERMAL MODEL 



The cut power law used in minijet fits and modified as a function of energy in the hadronization model pre- 
sented above accidentally coincides with the canonical Tsallis distribution. This distribution can be derived from a 
generalization of the Boltzmann entropy formula first suggested by Tsallis: 

= Y^T.iy^!-^^) (13) 

i 

In order to relate to our model we use q = 1 — 1/n. We get 

Sq = n^Wi{w~^^"' - 1) — > -^Wilnwi. (14) 

i i 

The canonical distribution maximizes the entropy with the following constraints: = 1 and WiSi = E. The 

use of a lagrange multipliers (5 leads to the following distribution which minimizes Sq — (3E: 

.,=ifi+^y" (15) 



Z \ n , 

with Z = + " for the proper normalization. 

Based on this distribution the free energy can be derived as well as relations reminding to rhe familiar thermody- 
namical relations. The Tsallis distribution can hence be regarded as a generalization of equilibrium thermodynamics. 
Among its physical realizations it is particularly interesting the case of finite heat reservoirs. 

In view of this formalism we interpret the one parton distribution as a canonical Tsallis distribution, 

wi = {l + j)-" = WTs.ms{E; T, 1 - 1/c) (16) 

with T = b/c. Recombination of n such partons leads to another Tsallis distribution with the same "temperature", 
but closer to the equilibrium case q = 1: 

Wn = Wi{E/nr = (1 + ^)-"" = WTs.ms{E; T, 1 - l/(nc)). (17) 



5 



Arriving at the finally observed hadrons with c « 8 and n = 2 . . . 6 one gets really close to the limiting exponential 

Woo = WTsams(^; T, 1) = e-^/^. (18) 

This property may explain the relative, and so far unexplained, success of the thermal model. 

Finally it is interesting to explore the potential of the Tsallis distribution with respect to particle numbers and the 
energy per particle, since these are the most profound predictions of the thermal model. For the sake of simplicity we 
consider here the analytically solvable case of massless particles. 

The number of a particle in a volume V with internal degeneracy factor d is given by the integral of the canonical 
Tsallis distribution. 

In the massless case, E = \p\, and we get 

In general this is a bigger number for finite c than for the Boltzmann distribution at nc ^ oo: N{c)/N{oo) > 1. This 
ratio is 2.4, 1.5, 1.3 for c = 8, 16, 24. 

The energy is given by another integral, 

which can as well be calculated in a closed form for massless particles, E = \p\: 

p _ ^ yiA ^ /oo'l 

27r2 (c- l)(c-2)(c-3)(c-4) ^ ' 

The energy per particle becomes E/N = ^ resulting in E/N w 1.07 GeV for b = 1.39 GeV and c = 7.9, and at the 
same time predicting a slope T w 176 MeV for the mr exponential at low momenta pt <^bcpa 10 GeV. 



V. CONCLUSION 



We reviewed limiting distributions of familiar short tailed and exceptional long tailed distributions. We have seen, 
that already a few, n = 3 — 4 independent random comopnents may lead to a distribution which is hard to distinguish 
from the limiting case on finite, fluctuating data samples. 

Non Gaussian limiting distributions exist with altered (non 1/y/n) scaling laws. For the high energy hadronization 
process it is particularly interesting to consider cut power law distributions, which in the limiting case combine to 
exponential spectra. This happens with the 1/n scaling, pointing out a possible interpretation as hadron spectra 
showing the Fourier transform of the distribution of the center of energy of the recombining partons. 

Minijet partons obtained using pQCD processes, and actually also experimentally observed spectra in pp collisions 
are Tsallis distributed at T w 170 — 200 MeV, and with E/N « 1 GeV. These two numbers can be simultanously 
described by a Tsallis distribution with the parameters b = 1.39 GeV and c = 7.9, when considering massless particles. 
The same simple consideration fails for the familiar Boltzmann distribution of the thermal equilibrium with a factor of 
2. At finite mass these integrals can be calculated only numerically. Considering massive particles with a mass aromid 
the pion mass, our results do not change dramatically. At the mass of 1 GeV, however, the Boltzmann distribution 
reaches E/N 1 GeV and the Tsallis distribution overshoots with a factor of 2. 
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